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Abstract 

We present a polymer quantization of the — A/r 2 potential on the positive real line and 
compute numerically the bound state eigenenergies in terms of the dimensionless coupling 
constant A. The singularity at the origin is handled in two ways: first, by regularizing 
the potential and adopting either symmetric or antisymmetric boundary conditions; second, 
by keeping the potential unregularized but allowing the singularity to be balanced by an 
antisymmetric boundary condition. The results are compared to the semiclassical limit of 
the polymer theory and to the conventional Schrodinger quantization on L2 0R+). The various 
quantization schemes are in excellent agreement for the highly excited states but differ for 
the low-lying states, and the polymer spectrum is bounded below even when the Schrodinger 
spectrum is not. We find as expected that for the antisymmetric boundary condition the 
regularization of the potential is redundant: the polymer quantum theory is well defined even 
with the unregularized potential and the regularization of the potential does not significantly 
affect the spectrum. 
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1 Introduction 



One of the most important outcomes expected of a successful theory of quantum gravity is a 
clear and unambiguous solution to the problems associated with the curvature singularities 
that are predicted by classical general relativity. This expectation is natural since quantum 
mechanics is known to cure classical singularities in other contexts, such as the hydrogen 
atom. 

In recent years there has been much work suggesting that Loop Quantum Gravity (LQG) 
[1] may indeed resolve gravitational singularities at least in the case of symmetry-reduced 
models, such as spatially homogeneous [5] and inhomogeneous [3] cosmologies and spherically 
symmetric black holes |H EJ |6]. Given the simplifications that these models entail, it is 
pertinent to ask which features of the LQG quantization scheme are crucial to the observed 
singularity resolution. 

There are two distinct, but related, features of the LQG quantization program that appear 
to play a role in achieving singularity resolution. The first is the fundamental discreteness that 
underlies LQG due to its focus on holonomies of connections and associated graphs embedded 
in a spatial manifold [7J . An analogous approach in a purely quantum mechanical context is 
so-called polymer quantization [3 19] , in which the Hamiltonian dynamics occurs on a discrete 
spatial lattice and the basic observables are the operators associated with location on the 
lattice and translation between lattice points. Polymer quantization provides a quantization 
scheme that is mathematically and physically distinct from Schrodinger quantization. 

The second apparently key ingredient in the LQG singularity resolution mechanism is the 
regularization of the singular terms in the Hamiltonian using a trick first introduced in this 
context by Thiemann [TJ. The regularization is achieved by first writing a classical inverse 
triad as the (singular) Poisson bracket of classical phase space functions whose quantum coun- 
terparts are known, and then defining the inverse triad operator as the commutator of these 
quantum counterparts. When applied to simple models this procedure gives rise to quantum 
operators with bounded spectra. The singularity is therefore kinematically "removed" from 
the spectra of relevant physical operators, such as the inverse scale factor. 

One question that arises concerns the role or perhaps the necessity of the Thiemann trick 
in singularity resolution in LQG. Recall that in the case of the hydrogen atom the singularity 
resolution is achieved by defining self-adjoint operators in a Hilbert space. This requires a 
careful choice of boundary conditions on the wave function [10] but does not require modi- 
fication of the singular 1/r potential. An example more relevant to quantum gravity is the 
reduced Schrodinger quantization of the "throat dynamics" of the Schwarzschild interior, 
which on imposition of suitable boundary conditions produced a discrete, bounded-from- 
below spectrum for the black hole mass 

Polymer quantization of the hydrogen atom was recently examined in [p2], retaining only 
the s-wave sector and regularizing the 1/r potential in a way that lets r take values on the 
whole real axis. The choice of symmetric versus antisymmetric boundary conditions at the 
singularity was found to have a signficant effect on the ground state even after the singularity 
itself had been regularized. In particular, in the limit of small lattice separation the ground 
state eigenenergy showed evidence of convergence towards the ground state energy of the 
conventionally-quantized Schrodinger theory only for the antisymmetric boundary condition. 
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In the present paper, we perform a similar polymer quantization of the more singular 1/r 2 
potential. When the potential is regularized, we shall find that the choice of the boundary 
condition again has a significant effect on the lowest-lying eigenenergies. However, we shall 
also find that the polymer theory with the antisymmetric boundary condition is well defined 
even without regularizing the potential, and with this boundary condition the regularized 
and unregularized potentials yield closely similar spectra. The boundary condition at the 
singularity is hence not only a central piece of input in polymer quantization, but it can even 
provide, along with the modification to the kinetic term, the pivotal singularity avoidance 
mechanism. While this is expected from general arguments that we will make explicit later 
on, it is interesting and reassuring to see the mechanism work in the special case of the 1/r 2 
potential, whose degree of divergence is just at the threshold where a conventional Schrodinger 
quantization will necessarily result into a spectrum that is unbounded below. The polymer 
treatment of this system is thus turning a Hamiltonian that is unbounded below into one 
that has a well-defined ground state. 

The 1/r 2 potential is interesting in its own right: it has a classical scale invariance that 
is broken by the quantum theory. In addition, this potential appears frequently in black 
hole physics, for example in the near horizon and near singularity behaviour of the quasinor- 
mal mode potential |13|. HH] , in the near horizon behaviour of scalar field propagation |15j 
and in the Hamiltonian constraint in Painleve-Gullstrand coordinates [16 . It may therefore 
conceivably be of direct relevance to quantum gravity. There is a substantial literature on 
Schrodinger quantization of this potential in L2(M+) (see for example [17 1 [18 ] [T9 l 120 1 [2T ] [22 ) |23] 
and the references in [24 J, but we are not aware of previous work on polymer quantization 
of this potential. 

Our paper is organized as follows. In Section [2] we review the Schrodinger quantization of 
the 1/r 2 porential in In Section [3] we formulate the polymer quantization of this 

system on a lattice of fixed size and describe the numerical method. We also include in this 
section a computation of the semi-classical polymer spectrum from the Bohr-Sommerfeld 
quantization condition, with a fixed polymerization length scale. The numerical results are 
presented in Section [4] and the conclusions are collected in Section [5j 



2 Schrodinger quantization 

We consider the classical Hamiltonian 



.2 A 



where the phase space is (r,p) with r > and A £ 1 is a constant. We shall take r, p and A all 
dimensionless, and on quantization we set h = 1. If physical dimensions are restored, r and 
p will become expressed in terms of a single dimensionful scale but A remains dimensionless. 
That the coupling constant is dimensionless is the speciality of a scale invariant potential. 



Quantization of H (2.1) is of course subject to the usual ambiguities. In particular, if 
one views H as an effective Hamiltonian that comes from a higher-dimensional configuration 
space via symmetry reduction, with r being a radial configuration variable, the appropriate 
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Hilbert space may be L2(M + ; m(r)dr), where m is a positive-valued weight function. If for 
example m(r) = r a , where a € M, then the ordering 

makes the quantum Hamiltonian H symmetric. If the wave function in L2(M + ; m{r)dr) is 
denoted by if), we may map ip to ij) € L2(M + ; dr) by ?/>(r) = r a / 2 ^(r), and i/ is then mapped 
in ^2(1^+; dr) to the Hamiltonian 

S = -—-±, (2.3) 

where 

X-A-?(!-l). (-> 

We shall consider any such mappings to have been done and take the quantum Hamiltonian 
to be simply 

_ pp. \ 
acting in the Hilbert space Z/2(M+; dr). 



To guarantee that the time evolution generated by H (2.5) is unitary, H must be specified 
as a self-adjoint operator on L2(M+; dr) [25]. A comprehensive analysis of how to do this was 
given in p2] (see also H3 ESI MH ED EH |23] ) . We shall review the results of [T9 in a way 
that displays the spectrum explicitly for all the qualitatively different ranges of A. 

Before proceeding, we mention that several recent quantizations of the 1 /r 2 potential first 
regularize the potential using various renormalization techniques |20| ED E2]. In particular, 
when the spectrum of a self-adjoint extension is unbounded below, these renormalization 
techniques need not lead to an equivalent quantum theory [23]. We shall here discuss only 
the self-adjoint extensions. 

To begin, recall [25] that the deficiency indices of H are found by considering normalizable 
solutions to the eigenvalue equation Hip = An elementary analysis shows that H is 

essentially self-adjoint for A < —3/4, but for A > —3/4 a boundary condition at r = 
is needed to make H self-adjoint. Physically, this boundary condition will ensure that no 
probability is flowing in/out at r = 0. 

2.1 A>l/4 

For A > 1/4, we write A = 1/4 + a 2 with a > 0. 

For E > 0, the linearly independent (non-normalizable) solutions to the eigenvalue equation 

Hip = Eip (2.6) 

are y/r J±i a (\/ n E r). These oscillate infinitely many times as r — > 0. To find the boundary 
condition, we consider the linear combinations 



il> E (r) ■= Vr" e^E~ ia / 2 J ia (VEr) + e ~ il3 E ia / 2 J- ia (VEr) , (2.7) 
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where ft is & parameter that a priori could depend on E. As ipE is periodic in j3 with 
period 2ir, and as replacing (3 by (3 + it multiplies ipE by —1, we may understand (3 periodic 
with period tt. For concreteness, we could choose for example (3 £ [0,7r). 
For the probability flux through r = to vanish, we need 



ipE d r ipE' — ipE' d r ipE — ► as r — > 



(2i 



for all E and E 1 , where the overline denotes the complex conjugate. Using the small argument 
asymptotic form (equation (9.1.7) in [26 ) 



(*/2)" 



as z 



0, 



(2.9) 



this is seen to require sin(/3 — (3') = 0, and hence (3 must be independent of E. The choice of 
the constant (3 hence specifies the boundary condition at the origin. 

To find the eigenvalues, we consider the normalizable solutions to (2.6). Such solutions exist 
only when E < 0, and they are y/r -Kj a (v 1 —E r). These solutions must satisfy at r — > the 
same boundary condition as ipE (2.7). Using the small argument asymptotic form (equations 
(9.6.2) and (9.6.7) in [26]) 



K V (Z) = K- V {z) 



7T 



(z/2Y 



(z/2Y 



r(-u + i) r(i/ + i) 



as z 



0, 



sin(^7r) 

this shows that the eigenvalues are 

E n = Eq exp(— 27rn/a), n G Z, 

where 

E = -exp[(2/3 + 7r)/a]. 

This spectrum is an infinite tower, with E n — > 0_ as n — > oo and £" n — > — oo as n 
The spectrum is unbounded from below. 

We note that Schrddinger quantization of a regulated form of the potential yields a semi- 
infinite tower of states that is similar to (2.11) as n — > oo but has a ground state |21j . The 
energy of the ground state goes to — oo when the regulator is removed. 

2.2 A = 1/4 



(2.10) 

(2.11) 

(2.12) 
* — oo. 



For A = 1/4, the solutions to the eigenvalue equation (2.6) for E > are y/rJo(vEr) and 
y/r Nq(\/E r) . We consider the linear combinations 



fe(r) := y/r { (cos/3) J (VEr) + (sin/3) 



TT 



No(VEr) - In 



£e 7 



Jo(VEr 



, (2.13) 



where 7 is Euler's constant and (3 is again a parameter that may be understood periodic 
with period tt and could a priori depend on E. As above, we find that (3 must be a constant 
independent of and its value determines the boundary condition at the origin. 
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Normalizable solutions to (2.6) exist only for E < 0. They are yfr Kq(\/— Er), and they 
must satisfy the same boundary condition as ipE (2.13) at r^O. Using the small argument 
expansion of Kq 



[26], we find that for (3 = there are no normalizable states, while for 
< (3 < 7r there is exactly one normalizable state, with the energy 



E = -4exp(-27 + 2cot/3). 



(2.14) 



2.3 -3/4 < A < 1/4 

For -3/4 < A < 1/4, we write A = 1/4 - v 2 with < v < 1 



The solutions to the eigenvalue equation (2.6) for E > are 
the linear combinations 



ijj E (r) ■= (cos P)E~ l//2 J l/ (VEr) + (sin f3)E u / 2 J- U {VE r 



J± v (yEr). Considering 



(2.15) 



we find as above that (3 is a constant, understood periodic with period ir, and its value 
specifies the boundary condition at the origin. 

Normalizable solutions to (2.6) exist only for E < 0. They are y/rK v {yf—Er) and must 
satisfy the same boundary condition as ipE (2.15) at r — > 0. Using the small argument 



asymptotic form of K v [SB], we find that there are no normalizable states for < j3 < vr/2, 
while for ir/2 < (3 < tt there is exactly one normalizable state, with the energy 



En 



-(-cot/3) 1/iy . 



We note that in the special case of a free particle, A 
trigonometric and exponential functions. 



(2.16) 

0, the Bessel functions reduce to 



2.4 A < -3/4 

For A < —3/4, we write A = 1/4 — v 2 with v > 1. H is now essentially self-adjoint. Any 



prospective normalizable solution to (2.6) would need to have E < and take the form 
\fr K u (\/—E r), but since now v > 1, these solutions are not normalizable and hence do not 
exist. 



3 Polymer quantization 

In this section we develop the polymer quantization of the 1/r 2 potential We proceed as 
in |12j . briefly reiterating the main steps for completeness. 

It is necessary to extend the r coordinate to negative values with the replacement r — > x G M 
in order to use central finite difference schemes at the origin. This will allow us to introduce 
at the origin both a symmetric boundary condition (with the regulated potential developed 



in subsection 3.2 ) and an antisymmetric boundary condition. 



The polymer Hilbert space on the full real line is spanned by the basis states 

<M*> = {w* (3 ' 1) 
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with the inner product 

(ipx,ipx') = 8 x ,x', (3-2) 
where the object on the right hand side is the Kronecker delta. The position operator acts 
by multiplication as 

(xijj)(x) = x-ip(x). (3-3) 

Defining a momentum operator takes more care. Consider the translation operator U, which 
acts as 

(U^)(x) = ^(x + fi). (3.4) 

In ordinary Schrodinger quantization we would have = e tfMp . Following [8], we hence define 
the momentum operator and its square as 

1 

2ijl ' 

4(2 -%-Ul), (3.5b) 
A* 



P = ^(Ufi-Ut), (3.5a) 



P 2 



where fi := (J./2. We may thus write the polymer Hamiltonian as 

Hp i = Tp i + V po \ , (3.6) 

where 

Tpol = -^(2-%-Ul), (3.7a) 
A* 

%oi = ~ (3.7b) 
Considering the action of x and Un, we see that the dynamics generated by (3.6) separates 



the polymer Hilbert space into an infinite number of superselection sectors, each having 
support on a regular /i-spaced lattice {A + nfi \ n £ Z}. The choice of {A | < A < fi} 
picks the sector. Since we wish to study singularity resolution, we concentrate on the A = 
sector, which we expect the singularity of the potenti al to affect most. We shall discuss the 
regularization of V po i at this singularity in subsection 



3.2 



3.1 Semiclassical polymer theory 

Before analyzing the full polymer quantum theory, we examine the semiclassical polymer 
spectrum using the Bohr-Sommerfeld quantization condition. 



Following [21 |5l E], we take the classical limit of the polymer Hamiltonian (3.6) by keeping 
the polymerization scale \i fixed and making the replacement Un — > e w , where p is the 
classical momentum. Note that this is different from the continuum limit in which fi goes to 
zero and the quantum theory is expected to be equivalent to Schrodinger quantization |27j . 

We assume A > 0. It follows, as will be verified below, that the classical polymer orbits 
never reach the origin, and we may hence assume the configuration variable x to be positive 
and revert to the symbol r. The classical polymer Hamiltonian takes thus the form 

_ sin 2 (^p) A 
#pol - — ^2 ^- (3-8) 
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Note that H po \ reduces to the classical non-polymerized Hamiltonian (2.1 ) in the limit p — > 0. 

A first observation is that the kinetic term in H po \ is non-negative and bounded above 
by 1/p 2 . Denoting the time-independent value of -ff po i on a classical solution by E, it follows 
that E is bounded above by 

E < £ max := 4> (3-9) 
and on a given classical solution r is bounded below by r > r_ , where 

--(AT- 

An elementary analysis shows every classical solution has a bounce at r = r_. For E > 
this is the only turning point, and the solution is a scattering solution, with r — ► 00 as 
i — * ±00. For E 1 < there is a second turning point at r = r + > r_, where 

and the solution is a bound solution, with r oscillating periodically between r+ and r_. Note 
that r+ is independent of p, and the outer turning point in fact coincides with the turning 
point of the non-polymerized classical theory. 

The classical polymer solutions are thus qualitatively similar to the classical non-polymerized 
solutions at large r, both for E > and for E < 0. What is different is that the polymer 
energy is bounded from above, and more importantly that the polymer solutions bounce at 
r = r_. In this sense the classical polymer theory has resolved the singularity at r = 0. The 
resolution depends on the polymerization scale: for fixed E, r_ = /2\/A[l + 0(fi 2 )] — » as 
p, — > 0, and for fixed /2, r_ — > py/X as E 1 — > 0. 

As the < solutions are periodic, we can use the Bohr-Sommerfeld quantization con- 
dition to estimate the semiclassical quantum spectrum. A subtlety here is that semiclassical 
estimates already in ordinary Schrodinger quantization with al/r 2 term involve a shift in the 
coefficient of this term [28 . Anticipating a similar shift here, we look at the Bohr-Sommerfeld 
quantization condition with A replaced by A e g, and we will then determine A e ff by comparison 
with the Schrodinger quantization. 



For a classical solution with given E, formula (3.8) implies (with A replaced by A e gf 



(3.12) 



\/sin 2 (pp) — p?E 
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Taking E < 0, the phase space integral J(E) := § r dp over a full cycle is hence 



J(E) = j> 



r dp 



2 V A e ff 



2 V A e ff 



7r/(2/Z) 

o \/sm 2 (p,p) — /2 2 ^ 

tt/2 



(p = y/li) 



dy 



o 



2\ 


/A e ff 




-/1 2 £ 


2\ 


/A c ff 




- p?E 



\/sm\y) - ^ 2 E 

dy 



1 - (1 - fi 2 E)~ 1 cos 2 (y) 



K[(l-fl 2 E) 



-1/2 



where K is the complete elliptic integral of the first kind 
expansion (8.113.3) in [29] yields 



(3.13) 

In the limit p?E — ► 0, the 



J(E) = 2VA 



cff 



In 



+ ' 



-E 



(3.14) 



^V-E, 

The Bohr-Sommerfeld quantization condition now states that the eigenenergies of the highly 



excited states are given asymptotically by J{E) = 2irn, where n 3> 1 is an integer. By (3.14) 
this gives the asymptotic eigenenergies 



En 



16 

^exp( 



n 



oo. 



(3.15) 



The Bohr-Sommerfeld estimate (|3.15|) agrees with the spectrum (|2.11|) obtained from con- 

A 



ventional Schrodinger quantization for A > 1/4, provided A e ff 
the self-adjoint extension for which 



and we choose in (2.11 ) 



P 



TT 



+ aln(4//x) (mod tt). 



(3.16) 



The shift A c ff = A — ^ is exactly that which arises in ordinary Schrodinger quantization of 
potentials that include a 1 jr 2 term: the reason there is the matching of the small r behaviour 
of the exact eigenstates to the WKB approximation. For a lucid analysis of this phenomenon 
in the quasinormal mode context, see the discussion between equations (23) and (28) in (28]. 
Note, however, that in our system the Bohr-Sommerfeld condition cannot be applied directly 
to the unpolymerized theory, since J(E) (3.13) diverges as p, — ► 0. 



3.2 Pull quantum polymer theory 



We now return to the full polymer quantum theory, with the Hamiltonian (3.6) and A G M. 

We write the basis states in Dirac notati on a s \mfj), where m S Z. Writing a state in this 
basis as ip = Ylm c m\ m ^) > it follows from (3.2) that the inner product reads = 
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^] m c m W Sm . The Hilbert space is thus ^(Z). It will be useful to decompose this Hilbert 
space as the direct sum Z/2(Z) = L^Z) © L^Z), where the states in the symmetric sector 
L|(Z) satisfy c m = c_ m and the states in the antisymmetric sector L^iZ) satisfy c m = — c_ m . 
The action of T po \ (3.7a) reads 



c m - c m+ i - c m _i) |m//} . 



(3.17) 



Tp i is clearly a bounded operator on L2(Z). T po \ is manifestly sy mmet ric, and an explicit 
solution of the eigenvalue equation T po \ip = Eip, given in equation ( 3.3l[ ) below, shows that 
there are no normalizable solutions with E = ±i. T po \ is hence essentially self-adjoint (|30j. 
Theorem X.2). It is also positive, since (^ip, T po \%l^j > for any ip ^ by the Cauchy-Schwarz 
inequality. 

The action of Vp \ (3.7b) reads 

^ P oi (^2 c m \rnfj,)j 

where 



(3.18) 



pol ._ 



(3.19) 



As (3.18) is ill-defined on any state for which c m ^ 0, V po \ is not a densely-defined operator 
on £2(2). We consider two ways to handle this singularity. 

The first way is to regulate V po [ explicitly. Recall that for x G M \ {0} we can write 



sgn(x) = 2 _(i(v / R) 



dx 

and on our lattice this can be implemented as the finite difference expression 

sgn(x) 1 



(3.20) 



(3.21) 



We hence define the regulated polymer version of sgn(x)/ \/\x\ by dropping the 0(fj, 2 ) term 
in (3.21), and we define the regulated polymer potential V*°f by raising this to the fourth 
power, 



A 



7 A 



or 



v4 



%m+l I 



J2fm s c m \mn) 



where 



(^m-v^ii) 4 . 



(3.22) 
(3.23) 
(3.24) 
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^pof ^ s c l ear ly a bounded essentially self-adjoint operator on L^ffi), and its operator norm is 
4|A|/(m 2 )- 

The regulated polymer Hamiltonian can now be defined by 



-"pol 



fpol + ^f. 



pol 



(3.25) 



It follows by the Kato-Rellich theorem ([30 , Theorem X.12) that H™^ is essentially self- 
adjoint on L2(Z) and bounded below by — 4|A|/(^ 2 ). Further, both T po \ and V™f leave 
and L?[(Z) invariant, and their boundedness and self-adjointness properties mentioned above 
hold also for their restrictions to anc ^ ^l(^)- ^ follows that H^j restricts to both 

L|(Z) and as a self-adjoint operator bounded below by — 4|A|/(// 2 ). We denote both 

of these restrictions by Hp^f, leaving the domain to be understood from the context. 



The second way to handle the singularity of V po i (3.18) is to restrict at the outset to the 
antisymmetric subspace L^Z), where V po \ is essentially self-adjoint and its operator norm is 
|A|/(/x 2 ). It follows as above that the unregulated polymer Hamiltonian 

Hpoi = r po i + V P oi 

on ^(Z) is essentially self-adjoint and bounded below by — |A|/(// 2 ). 

Two comments are in order. First, H veg can be written in terms of operators as 



(3.26) 



h 1 H 1 V 



The potential in (3.27) can hence be viewed as arising by the substitution 

sgn(x) 



2 



ul{VW\A} 



(3.27) 



(3.28) 



in place of (3.20). This method is similar to Thiemann's regularization of inverse triad 



operators in loop quantum gravity [7] . 

Second, the regulated potential vanishes at the origin but is greater in absolute value than 
the unregulated potential for \m\ > 1. However, the difference is significant only for the 
lowest few \m\, and the two potentials quickly converge as \m\ — > oo, as shown in Figure [TJ 
The regulated and unregulated potentials hence differ significantly only near the singularity. 



3.3 Eigenstates and the numerical method 

We are now ready to look for the eigenstates of the Hamiltonian. Writing the eigenstate 
as V = Em*™ \mfj,) and denoting the eigenvalue by E, the regulated eigenvalue equation 
H^tp = Eip an d the unregulated eigenvalue equation H po \ip = Etp both give a recursion 
relation that takes the form 

Cm (2 - fi 2 E - Xf m ) = c m+ i + Cm-!, (3.29) 
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where f m = fm g (3.24) for the regulated potential and f m = fm 1 (3.19) for the unregulated 
potential. Note that the polymerization scale \x enters this recursion relation only in the 
combination [i 2 E, whether or not the potential is regulated. This is a direct consequence of 
the scale invariance of the potential. 

From now on, we take A > and E < 0. 

We use the "shooting method" that was applied in [12] to the polymerized l/\x\ potential. 



For large |m|, (3.29) is approximated by 



Cm (2 - Li 2 E) = c m+ i + Cm-i. 



(3.30) 



The linearly independent solutions to ( 3.30[ ) are 




^ 2 



±m 



(3.31) 



The upper (respectively lower) sign gives coefficients that increase (decrease) exponentially 
asm-> oo. We can therefore use (3.31) with the lower sign to set the initial conditions at 
large positive m [51] . 

To set up the shooting problem, we choose a value fo r n 2 E and begin with some uiq 3> 
A to find 



and 



c, 



mo-l 



using the approximation (3.31). We then iterate downwards 



with (3.29). In the antisymmetric sector, we stop the iteration at cq and shoot for values 
of [i 2 E for which cq = 0. This shooting problem is well defined both for the unregulated 



potential (3.18) and for the regulated potential (3.23), since the computation of cq via (3.29) 
does not require evaluation of f m at m = 0. In the symmetric sector, we stop at the iteration 
at e_i and shoot for values of fJ?E for which c_i = c\. As the computation of c_i requires 
evaluation of f m at m = 0, the symmetric sector is well defined only for the regulated 
potential. 



4 Results 

We shall now compare the spectra of full polymer quantization, Bohr-Sommerfeld polymer 
quantization and ordinary Schrodinger quantization. We are particularly interested in the 
sensitivity of the results to the choice of the symmetric versus the antisymmetric sector. 

First of all, we find that when the potential is regulated, the choice of the symmetric versus 
antisymmetric boundary condition in the full polymer quantum theory has no significant 
qualitative effect for sufficiently large A, the only difference being slightly lower eigenvalues 
for the symmetric boundary condition. The lowest five eigenvalues in the two sectors are 
shown in Table [T] for A = 2. This is in a sharp contrast with what was found in [pj] for the 
1 jr potential, where the symmetric sector contained a low-lying eigenvalue that appeared to 
tend to — oo as the polymerization scale was decreased. 

Another key feature is that for sufficiently large A there is indeed a negative energy ground 
state. For 3 < A < 4, the plots of the lowest eigenvalues as a function of A in Figures [2] and 
[3] show that the analytic lower bound obtained in subsection 3.2 is accurate within a factor 
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antisymmetric 


symmetric 


E 


-6.14 


-6.37 


Ei 


-2.35 • 10~ 2 


-2.43- 1(T 2 


E 2 


-2.03 • 10~ 4 


-2.10- 10" 4 


E 3 


-1.76 • 1(T 6 


-1.82 • 10~ 6 


E± 


-1.52 • 10" 8 


-1.57- 10" 8 



Table 1 : The lowest five eigenvalues of the regulated potential with antisymmetric and sym- 
metric boundary conditions (A = 2, fj, = 1). 



of 1.2 for the regulated potential in the antisymmetric sector and within a factor of two for 
the unregulated potential. 

Figures [2] and [3] also indicate that the lowest eigenvalues converge towards zero as A de- 
creases, for both the unregulated and regulated potentials, with the unregulated eigenvalues 
reaching zero slightly before the regulated. Near E n = the relationship is quadratic in A 
while the plots straighten out to a linear relationship for larger A. 

The numerics become slow as the energies are close to zero. We were unable to investigate 
systematically whether bound states exist for A < 1 /4, and in particular to make a comparison 
with the single bound state that occurs in Schrodinger quantization with certain self-adjoint 
extensions. For A slightly below 1/4, we do find one bound state, but we do not know 
whether the absence of further bound states is a genuine property of the system of an artefact 
of insufficient computational power. This would be worthy of further investigation. The 
eigenvalues show a similar dependence on A for both regulated and unregulated potentials, 
with the energies for the regulated potential being lower than those for the unregulated 
version as one would expect from comparing the potentials as in Figure [TJ 

For A > 1/4, we find that the eigenvalues E n depend on n exponentially, except for the 
lowest few eigenvalues (n = 0, 1). The coefficient in the exponent is in close agreement with 
the exact Schrodinger spectrum (2.11) and with the Bohr-Sommerfeld polymer spectrum 
( 3.15[ ) with A c ff = A — 1/4. Representative spectra are shown as semi-log plots in Figures [4] 
and [5] where the linear fit is computed using only the points with n > 2. By matching the 
linear fit to the Schrodinger spectrum (2.11 ) and reading off the self-adjointness parameter /?, 
we can determine the self-adjoint extension of the Schrodinger Hamiltonian that matches the 
polymer theory for the highly-excited states. The results, shown in Figures [6] and [7J show 
that the self-adjointness parameter (3 depends linearly on the coupling parameter a, and the 
slope in this relation is within 10 per cent of the slope obtained from the Bohr-Sommerfeld 
estimate ( |3 16| ), In 8 » 2.0794 (for /j, = l,p, = 1/2). 

Finally, our numerical eigenvalues E n are in excellent agreement with the analytic approxi- 
mation scheme of [16] , provided this scheme is understood as the limit of large A with fixed n. 
If the numerical results shown in figures [6] and [7] are indicative of the complementary limit 
of large n with fixed A, they show that the approximation scheme of [16] does not extend to 
this limit. 
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5 Conclusions 



We have compared Schrodinger and polymer quantizations of the 1/r 2 potential on the pos- 
itive real line. The broad conclusion is that these quantization schemes are in excellent 
agreement for the highly-excited states and differ significantly only for the low-lying states. 
In particular, the polymer spectrum is bounded below, whereas the Schrodinger spectrum 
is known to be unbounded below when the coefficient of the potential term is sufficiently 
negative. We also find that the Bohr-Sommerfeld semi classical quantization condition repro- 
duces correctly the distribution of the highly-excited polymer eigenvalues. At some level this 
agreement is not surprising, since one expects that for any mathematically consistent quan- 
tization scheme, in some appropriate large n, semi-classical limit, the spectra should agree. 
For anti-symmetric boundary conditions both Schrodinger and regulated and unregulated 
polymer obey the criteria, so it is perhaps not surprising that they agree at least for energies 
close to zero. It is somewhat surprising that they agree so well for low n (where "low" is in 
the context of the polymer spectra which are bounded below) . 

A central conceptual point was the regularization of the classical r = singularity in the 
polymer theory. We did this first by explicitly regulating the potential, using a finite differenc- 
ing scheme that mimics the Thiemann trick used with the inverse triad operators in LQG [7|: 
this method allows a choice of either symmetric or antisymmetric boundary conditions at the 
origin. We then observed, as prevously noted in |16j . that the singularity can alternatively be 
handled by leaving the potential unchanged but just imposing the antisymmetric boundary 
condition at the origin. The numerics showed that all of these three options gave very similar 
spectra, and the agreement was excellent for the highly-excited states. 

To what extent is the agreement of these three regularization options specific to the 1/r 2 
potential? Consider the polymer quantization of the Coulomb potential, —1/r. When the 
Coulomb potential is explicitly regulated, it was shown in [12] that the choice between the 
symmetric and antisymmetric boundary conditions makes a significant difference for the 
ground state energy. We have now computed numerically the lowest five eigenenergies for 
the unregulated —1/r potential with the antisymmetric boundary condition, with the results 
shown in Table [2j Comparison with the results in (12] shows that the regularization of the 
potential makes no significant difference with the antisymmetric boundary condition. As 
noted in [12 , for sufficiently small lattice spacing the antisymmetric boundary condition 
spectrum tends to that which is obtained in Schrodinger quantization with the conventional 
hydrogen s-wave boundary condition (10J . 

We conclude that in polymer quantization of certain singular potentials, a suitably-chosen 
boundary condition suffices to produce a well-defined and arguably physically acceptable 
quantum theory, without the need to explicitly modify the classical potential near its singu- 
larity: the antisymmetric boundary condition effectively removes the r = eigenstate from 
the domain of the operator 1/r 2 by requiring cq = in the basis state expansion Y^m c m\mij)- 
A similar observation has been made previously in polymer quantization of a class of cosmo- 
logical models, as a way to obtain singularity avoidance without recourse to the Thiemann 
trick [32[ [33] . and related discussion of the self-adjointness of polymer Hamiltonians arising 
in the cosmological context has been given in |34] . While we are not aware of a way to relate 
our system, with the —1/r 2 potential and no Hamiltonian constraints, directly to a specific 
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n 


E n 





-0.250 


1 


-0.0625 


2 


-0.0278 


3 


-0.157 


4 


-0.0100 



Table 2: The lowest five eigenvalues for the unregulated Coulomb potential with the anti- 
symmetric boundary condition (A = 1, /j, = 0.1). 



cosmological model, it is nonetheless reassuring that the various techniques we have used in 
this case for dealing with the singularity all lead to quantitatively similar spectra. Whether 
this continues to hold in polymer quantization of theories that are more closely related to 
LQG is an important open question that is currently under investigation [35 . 



Acknowledgements 

We thank Jack Gegenberg, Sven Gnutzmann and Viqar Husain for helpful discussions and 
correspondence. GK and JZ were supported in part by the Natural Sciences and Engineering 
Research Council of Canada. JL was supported in part by STFC (UK) grant PP/D507358/1. 



References 



[1] T. Thiemann, Modern Canonical Quantum General Relativity (Cambridge University 
Press, Cambridge, 2007). 

[2] A. Ashtekar, T. Pawlowski and P. Singh, "Quantum nature of the big bang: Im- 
proved dynamics," Phys. Rev. D 74, 084003 (2006) [arXiv:gr-qc/0607039] ; A. Ashtekar, 
T. Pawlowski, P. Singh and K. Vandersloot, "Loop quantum cosmology of k = 1 FRW 
models," Phys. Rev. D 75, 024035 (2007) |arXiv:gr-qc/0612104| ; K. Vandersloot, "Loop 
quantum cosmology and the k = -1 RW model," Phys. Rev. D 75, 023523 (2007) 
[arXiv:gr- qc/0612070 ; M. V. Battisti, O. M. Lecian and G. Montani, "Polymer quan- 
tum dynamics of the Taub universe," Phys. Rev. D 78, 103514 (2008) [arXiv:0806.0768| 
[gr-qc]]. 

[3] M. Martin-Benito, L. J. Garay and G. A. Mena Marugan, "Hybrid Quantum Gowdy 
Cosmology: Combining Loop and Fock Quantizations," Phys. Rev. D 78, 083516 (2008) 



|arXiv:0804. 1098 [gr-qc]]. 



[4] V. Husain and O. Winkler, "Quantum Hamiltonian for gravitational collapse," Phys. 
Rev. D 73, 124007 (2006) |arXiv:gr-qc/0601082] . 

[5] A. Ashtekar and M. Bojowald, "Quantum geometry and the Schwarzschild singularity," 
Class. Quant. Grav. 23, 391 (2006) [arXiv:gr-qc/0509075l . 



14 



[6] L. Modesto, "Loop quantum black hole," Class. Quant. Grav. 23, 5587 (2006) |arXiv:gr- 



qc/0509078 ; "Black hole interior from loop quantum gravity," |arXiv:gr-qc / 061 1043 



C. G. Boehmer and K. Vandersloot, "Loop Quantum Dynamics of the Schwarzschild 
Interior," Phys. Rev. D 76, 104030 (2007) |arXiv:0709.2129 [gr-qc]]; "Stability of the 
Schwarzschild Interior in Loop Quantum Gravity," Phys. Rev. D 78, 067501 (2008) 
[arXiv:0807.3 042 [gr-qc]]; M. Campiglia, R. Gambini and J. Pullin, "Loop quantization of 
spherically symmetric midi-superspaces: the interior problem," AIP Conf. Proc. 977, 52 
(2008) |arXiv:0712.0817 [gr-qc]]; R. Gambini and J. Pullin, "Black holes in loop quantum 
gravity: the complete space-time," Phys. Rev. Lett. 101, 161301 (2008) [arXiv:0805.1187| 
[gr-qc]]. 

T. Thiemann, "Quantum spin dynamics (QSD)," Class. Quant. Grav. 15, 839 (1998) 
|arXiv:gr-qc/9 606089 . 

A. Ashtekar, S. Fairhurst and J. L. Willis, "Quantum gravity, shadow states, and quan- 
tum mechanics," Class. Quant. Grav. 20, 1031 (2003) [arXiv:gr-qc/0207106l . 

H. Halvorson, "Complementarity of representations in quantum mechanics," Studies 
Hist. Philos. Mod. Phys. 35, 45 (2004) [arXiv:quant-ph/0110102"] . 

C. J. Fewster, "On the energy levels of the hydrogen atom," |arXiv:hep-th /9305102, 

J. Louko and J. Makela, "Area spectrum of the Schwarzschild black hole," Phys. Rev. 
D 54, 4982 (1996) [arXiv:gr-qc/9605058| . 

V. Husain, J. Louko and O. Winkler, "Quantum gravity and the Coulomb potential," 
Phys. Rev. D 76, 084002 (2007) |arXiv:0707.0273 [gr-qc]]. 

L. Motl, "An analytical computation of asymptotic Schwarzschild quasinormal frequen- 
cies," Adv. Theor. Math. Phys. 6, 1135 (2003) [arXiv:gr-qc/0212096] . 

G. Kunstatter, "d-dimensional black hole entropy spectrum from quasi-normal modes," 
Phys. Rev. Lett. 90, 161301 (2003) [arXiv:gr-qc/0212014| . 

H. E. Camblong and C. R. Ordonez, "Black hole thermodynamics from near-horizon 
conformal quantum mechanics," Phys. Rev. D 71, 104029 (2005) |arXiv:hep-th/041 1008 . 

J. Gegenberg, G. Kunstatter and R. D. Small, "Quantum structure of space near a black 
hole horizon," Class. Quant. Grav. 23, 6087 (2006) |arXiv:gr-qc/0606002| . 

K. M. Case, "Singular Potentials," Phys. Rev. 80, 797 (1950). 

W. M. Frank, D. J. Land and R. M. Spector, "Singular Potentials," Rev. Mod. Phys. 
43, 36 (1971). 

[19] H. Narnhofer, "Quantum theory for l/r 2 -potentials," Acta Phys. Austriaca 40, 306 
(1974). 



15 



[20] K. S. Gupta and S. G. Rajeev, "Renormalization in quantum mechanics," Phys. Rev. D 
48, 5940 (1993) |arXiv:hep-th/9305052| . 

[21] H. E. Camblong, L. N. Epele, H. Fanchiotti and C. A. Garcia Canal, "Renormalization of 
the inverse square potential," Phys. Rev. Lett. 85, 1590 (2000) [arXiv:hep-th/0003014] ; 
"Dimensional transmutation and dimensional regularization in quantum mechanics. I: 
General theory," Annals Phys. 287, 14 (2001) [arXiv:hep-th/0003255| . 

[22] H. E. Camblong and C. R. Ordonez, "Renormalization in conformal quantum mechan- 
ics," Phys. Lett. A 345, 22 (2005) |arXiv:hep-th/0305035l ; "Anomaly in conformal quan- 
tum mechanics: From molecular physics to black holes," Phys. Rev. D 68, 125013 (2003) 
[arXiv:hep-th/0303166] ; G. N. J. Ananos, H. E. Camblong and C. R. Ordonez, "SO(2,l) 
conformal anomaly: Beyond contact interactions," Phys. Rev. D 68, 025006 (2003) 
[arXiv:hep-th/0302197] . 

[23] H. E. Camblong, L. N. Epele, H. Fanchiotti, C. A. Garcia Canal and C. R. Ordonez, "On 
the Inequivalence of Renormalization and Self- Adjoint Extensions for Quantum Singular 
Interactions," Phys. Lett. A 364, 458 (2007) | arXiv:hep-th/0604018] . 

[24] T. Fulop, Symmetry, Integrability and Geometry, Proceedings of the 3rd 
Microconference "Analytic and Algebraic Methods III", 
http:/ / www.emis.de/journals / SIGMA /Prague2007.html 

[25] W. E. Thirring, Quantum mathematical physics, 2nd edition (Springer, Berlin, 2002). 

[26] Handbook of Mathematical Functions, edited by M. Abramowitz and I. A. Stegun (Dover, 
New York, 1965). 

[27] A. Corichi, T. Vukasinac and J. A. Zapata, "Polymer Quantum Mechanics and its Con- 
tinuum Limit," Phys. Rev. D 76, 044016 (2007) [arXiv:0704.0007| [gr-qc]]. 

[28] N. Andersson and C. J. Howls, "The asymptotic quasinormal mode spectrum of non- 
rotating black holes," Class. Quant. Grav. 21, 1623 (2004) [arXiv:gr-qc/0307020] . 

[29] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products, 4th edition 
(Academic, New York, 1980). 

[30] M. Reed and B. Simon, Methods of Modern Mathematical Physics II: Fourier Analysis, 
Self-adjointness (Academic, New York, 1975). 

[31] S. Elaydi, An Introduction to Difference Equations , 3rd edition (Springer, New York, 
2005). 

[32] A. Ashtekar, A. Corichi and P. Singh, "On the robustness of key features of loop quantum 
cosmology," Phys. Rev. D 77, 024046 (2008) [arXiv:0710.3565| [gr-qc]]. 

[33] A. Corichi and P. Singh, "Is loop quantization in cosmology unique?," Phys. Rev. D 78, 
024034 (2008) |arXiv-0805.0136 [gr-qc]]. 



16 



[34] W. Kaminski and J. Lewandowski, "The flat FRW model in LQC: the self-adjointness," 
Class. Quant. Grav. 25, 035001 (2008) [arXiv:0709.3l20 [gr-qc]]. 

[35] G. Kunstatter, J. Louko and J. Ziprick, in progress. 



17 




X 



Figure 1: The solid (red) line is the regulated — 1/x 2 potential (A = 1, /J, = 1). The dashed 
(blue) line is the unregulated potential. 
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Figure 2: The lowest two energy levels as a function of A for the unregulated potential with 
antisymmetric boundary conditions. 
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Figure 3: The lowest two energy levels as a function of A for the regulated potential with 
antisymmetric boundary conditions. 
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Figure 4: ln(E n ) vs. n for the unregulated potential with antisymmetric boundary conditions 
with linear fits (R 2 = 1). 




Figure 5: ln(E n ) vs. n for the regulated potential with antisymmetric boundary conditions 
with linear fits (R 2 = 1). 
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Figure 6: [3 vs. a with /j, = 1 for the unregulated potential with antisymmetric boundary 
conditions with a linear fit (R 2 = 1). 
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Figure 7: (3 vs. a with \i = 1 for the regulated potential with antisymmetric boundary 
conditions with a linear fit (R 2 = 1). 
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